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Abstract 

Investigating the Z3 symmetry in Quantum Chromodynamics (QCD) we 
show that fuh QCD with a vacuum of vanishing baryonic number does not 
lead to metastable phases. Rather in QCD with dynamical fermions, the 
degeneracy of Z3 phases manifests itself in observables without open triality. 
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The Lagrangian in pure gluonic Quantum Chromodynamics (QCD) has global Z3 sym- 
metry. It is very well known that in the low temperature phase Z3 symmetry leads to 
confinement of color charges and in the high temperature phase Z3 symmetry is sponta- 
neously broken, color charges are screened by the gluon field. In QCD with dynamical 
quarks the fermionic contribution to the Lagrangian breaks Z3 symmetry explicitly. 

Recently, in refs. it has been concluded on the one hand that this explicit Z3 symme- 
try breaking leads to the existence of metastable states at arbitrary high temperatures. On 
the other hand the authors of refs. believe that Zn phases cannot be prepared as real 
macroscopic systems. In the following we will try to remove some confusion and show that 
QCD with dynamical fermions has degenerate Z3 phases and therefore also ordered-ordered 
phase transitions like pure gluonic QCD. Furthermore, we emphasize that QCD is able to 
describe symmetric systems only, i.e. systems of zero triality. In order to introduce 
our concept we study first the pure gluonic system and then turn to QCD with dynamical 
fermions. 

In lattice QCD the gluon field appears in the form of SU{3) matrices Ux,^ which are 
defined on links (x, /i) of a four-dimensional euclidean lattice. In finite temperature lattice 
QCD the correlation function < L{fi) ■ ■ ■ L*{f]\f) > of several Polyakov loops 

1 Nt 

m = TT Tr n U(F,t),o (1) 
t=l 

can be connected 0] with the free energy F of N infinitely heavy quarks q or antiquarks q 
at the corresponding positions ri, • ■ ■ , tat at temperature T relative to the free energy of the 
vacuum 

F{q{n), ■ ■ ■ , g(r-V)) = -Tin < L(ri) ■ ■ ■L*(rN) > . (2) 
The thermodynamical average < • ■ • > is computed using the partition function 

Z = jv[U] e-^[^l (3) 

as a "sum" over all gauge field configurations U. In pure gluonic QCD for most problems 
the appropriate choice for S is the Wilson action 
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Sg[U] =(3Y^ fl - Tr U,,,^ , = 

where the plaquettes Ux,^v are built of four hnks in the /xz^-plane of a four dimensional 
euclidean lattice. 

The action <S'G[t^] is symmetric. This means that multiplication of all links in direction 
/i = in the three dimensional x,y,z-torus with fixed t, e.g. t = by a Z3 element leaves 
the action invariant. A single Polyakov loop L(r) transformes under non-trivially (its 
triality T is one), and therefore the distribution of L values can be used as an indicator for 
spontaneous breaking of Z^ symmetry. 

Let us now discuss the results of lattice calculations and their common interpretation. 
Monte-Carlo calculations in lattice QCD show usually a characteristic behaviour of the 
spatial average L of L(r) during a Monte-Carlo simulation. In the confinement regime 
(low 13) L scatters symmetrically around zero in the complex plane. In the deconfined 
phase there appear three Z^ symmetric maxima of the Polyakov loop distribution in 0°- and 
±120°-directions. As the number of tunneling events between the maxima decreases with 
increasing (3 one commonly expects L to be frozen in the thermodynamical limit in one of 
the Z3 directions and thus to get spontaneous Z3 symmetry breaking. The appearance of 
the three peak structure of the L-distribution in the deconfinement regime is considered as 
a demonstration that spontaneous symmetry breaking on a finite lattice can never happen 
exactly. In the thermodynamical limit one may obtain arg < L >= 1,±^. Besides the 
unpleasant fact that for arg(< L >) = ±^ one should reconsider the definition (0) the 
above interpretration of Monte-Carlo results leads in full QCD to doubtful consequences. 
Therefore, we prefer another point of view: As / T>\U] e"'^'''^^ is Z3 invariant and the 
Polyakov loop (|I|) is not, one has to get 

< L >= (5) 

in the confined and also in the deconfined phase. There is no contradiction between this 
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statement and the spontaneous breaking of Z3 symmetry in the deconfined phase. It rather 
shows that the expectation value < L > is not the appropriate observable for the demon- 
stration of spontaneous breaking of Z3 symmetry. One should rather use < L{0)L*{r) > 
correlations or the (Z3 symmetric) L distributions. The latter is equivalent to the determi- 
nation of Polyakov loops in Z3 invariant representations (octet, decuplet, antidecuplet, etc). 
After this rather trivial observation that for a Z3 invariant Lagrangian observables of non 
vanishing triality have expectation value zero we will turn to the more interesting case of 
full QCD. 

In full QCD with dynamical fermions the action S contains a fermionic contribution Sp 

S = Sq -\- Sp, 

x,x' 

which breaks symmetry explicitly. In the Kogut-Susskind formulation the fermionic 
matix M reads 



' r) ^ y X,y,Ux,^^x' ,X+^ ^ x' ,llU^/ ^6x' ^x—fl) ~l~ ^^X,x' ; . 

(7) 



P _ (_1\X1+X2^ \-Xfj,-l 

^ x,fi — \ ^) ; 



with one-component Grassmann variables ip^ and ip^- The integration over the Grassmann 
variables in the path integral 

Z = jv[U,ij,i;]e-' (8) 

can be performed analytically and leads to the fermionic determinant detM 

Z = jv[U] e-'^^t^'detM. (9) 

It is easy to see that the fermionic action violates Z-^ symmetry 0. In the usual MC- 
iterations one finds that in the low temperature phase the expectation value < L > of a 
Polyakov loop L is a small positive number and in the high temperature phase the three 
peak structure of the distribution becomes asymmetric. The maximum at 0° is favoured 
and the two lower maxima at angles of around ±120° are symmetric. This fact is usually 
interpreted as a consequence of the violation of Z^ symmetry. 



Effective models based on tliis Z3 asymmetric L-distribution are discussed in refs. HJ-Q. 
As mentioned in tlie introduction tliey predict metastable states, but witli unpliysical prop- 
erties. 

In order to resolve this contradiction we start from the statement that the baryonic num- 
ber B of the QCD vacuum should be zero. It follows that also the triality T is zero. Instead 
of the canonical ensemble with B = T = which should be well suited, a grand canonical 
description is allowed if it predicts the same behaviour in the thermodynamical limit. Com- 
mon lattice Monte-Carlo calculations in full QCD use the full fermionic determinant detM 
and thus simulate a grand canonical ensemble with chemical potential = containing all 
three triality sectors. There appears now the following loophole in the usual argumentation. 
The determination of the Polyakov loop expectation value in full QCD gives < L and 
a Z3 asymmetric L-distribution. In assuming to have measured the L-distribution of the 
vacuum several authors use effective potentials for the time component of the gluon field 
arising from one loop effects and predict metastable states |]I],|2| or find properties of these 
states 1^,^ which are not possible in Minkowski space. However, it has been overlooked that 
the determinant detM may be decomposed in three different triality sectors. The decom- 
position can easily be derived in the following way: Via a hopping parameter expansion the 
fermionic determinant can be decomposed into contributions according to the number of 
closed loops of gauge field U which wind around the 4-dimensional lattice in time direction 
||8HlO|. We can then define the triality operator T (with eigenvalues 0, ±1) which determines 
the triality of such closed loops, and with his help a projection operator Pj- on triality T 



3 q=0,±l 

Then the triality zero contribution of the fermionic determinant is 



(10) 



^ q=0,±l 

and more general a projection operator on triality T f\ 




(11) 



§The implementation of this projection operator in the lattice formulation is, at least in principle, 
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,q2m{f-T)/'i 



(12) 



and in general 



detM = detrM, 



(13) 



r=o,±i 



with 



detrM = PrdetM. 



(14) 



From eq. ([131) it is obvious that the observable L having triality T = 1 tests only the 
T = — 1 sector of the fermionic determinant which is not the vacuum sector. Using the 
correct T = vacuum sector for the fermionic determinant results in a symmetric L- 
distribution, an expectation value < L >= and no prediction of metastable and unphysical 
vacuum states. Let us discuss this in more detail. As the measure 'D\U] and the gluonic 
Lagrangian S'^ft/] are Z3 symmetric, the product detM[?7]0r[^] in 



has to be also symmetric to get < Or >7^ 0. In other words Z3 violating contributions 
in detM [?7]07-[[/] are automatically eliminated by the path integral. Only the triality — T 
component det_7-M of the fermionic determinant survives. det_rM has T quark loops less 
than antiquark loops winding around the lattice in time direction. 

Applying this statement to the example of the single Polyakov loop L representing an 
infinitely heavy quark Q, the fermionic determinant supplies a triality —1 state of light 
quarks which will mostly consist of a light antiquark q in order to color neutralize the heavy 
quark Q. A finite value of < L > in full QCD therefore does not mean that Z^ symmetry is 



very simple. One has to multiply all links in time direction for an arbitrary time step with the 
phase e''^'^*/^ and to determine with this new link variables the determinant. The result has to be 
multiplied with e'^'^'^^'^ 1'^ . Finally, the sum over phases g = 0, ±1 has to be executed. 




(15) 
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broken but rather that the fermionic system is flexible enough to neutralize external charges. 
The Polyakov expectation value < L > is a thermodynamical mixture 

<L>= ... + e-^(Q^")/^ + e-^(««^)/^ + ... (16) 

of a heavy-light meson Qq, a heavy-light baryon Qqq, ... . This heavy-light states should be 
bound in the hadronic and "ionized" in the quark- gluon plasma phase. The charge density 
of such light quarks around static sources has been measured in ref. [|ll[ . 

For the system of a heavy quark Q at position ri and a heavy antiquark Q at position r2 

Lin)L%f2) (17) 

a totaly different component of the fermionic determinant, the triality zero component, 
contributes. It is evident that this component is Z-^ symmetric, i.e. the three transformed 
sectors which differ in L by factors contribute with the same Boltzmann factor. This 
example shows that for observables of triality zero full QCD is Z^ symmetric. The fermionic 
action does not destroy Z^ symmetry as is the general belief. 

We can even extend this statement for triality non zero observables. The above discussion 
for a single Polyakov line has demonstrated that for Z3 violating observables full QCD 
restores the Z3 symmetry by compensating them to triality zero. Therefore, the effective 
models of refs. are effective models of heavy-light hadrons in the QCD vacuum and 

not effective models of "the" QCD vacuum. 

We would like to summarize that the vacuum of QCD with dynamical fermions has 
triality zero and therefore degenerate Z3 phases and ordered-ordered phase transitions like 
pure gluonic QCD. The Polyakov loop can still be used as an order parameter for the breaking 



of Z3 symmetry. There is no need to search new order parameters |]T2| for full unquenched 
QCD. In usual lattice Monte-Carlo simulations this triality zero sector can be projected 
out by using triality zero obsevables only. Another method would be to use detoM as the 
fermionic Boltzmann factor. However, an implementation of a projection to fixed triality T 
(e.g. T = 0) will be rather complicated due to the fact that usually not detM is simulated 
but (detM)^ whose T = component 



(detM)^^o = (detoM)2 + 2detiMdet_iM (18) 

contains contributions from triality ±1 components. Therefore, one will have to rely on 
common lattice calculations with the limitations that one should not be mislead by the 
mixing of different triality sectors in the grand canonical ensemble. 
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After finishing the first version of this article we noticed that J. Polonyi [|l3l introduced 



a concept similar to ours. The main conceptual difference to Polonyi's idea is that we 
demonstrate that common Monte- Carlo calculations use the correct Boltzmann factor for 
the fermionic vacuum if only triality zero observables are measured. 
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